Potential inversion with subbarrier fusion data revisited 
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We invert experimental data for heavy-ion fusion reactions at energies well below the Coulomb 
barrier in order to directly determine the internucleus potential between the colliding nuclei. In con- 
trast to the previous applications of the inversion formula, we explicitly take into account the effect 
of channel couplings on fusion reactions, by assuming that fusion cross sections at deep subbarrier 
energies are governed by the lowest barrier in the barrier distribution. We apply this procedure 
to the 16 + 144 Sm and le O + 208 Pb reactions, and find that the inverted internucleus potential 
are much thicker than phenomenological potentials. A relation to the steep fall-off phenomenon of 
fusion cross sections recently found at deep subbarrier energies is also discussed. 
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Nuclear reactions are primarily governed by the 
nucleus-nucleus potential. In particular, the Coulomb 
barrier, which appears due to a strong cancellation be- 
tween the attractive nuclear force and the long-range 
repulsive Coulomb interaction, plays a decisive role in 
heavy-ion collisions. Several methods have been pro- 
posed to compute the real part of the internuclear po- 
tential. Among them, the double folding model has 
been often employed and has enjoyed a success in de- 
scribing elastic and inelastic scattering for many systems 
0, S El- The Woods-Saxon form has also often been 
used to parametrize the inter-nuclear potential [4]. The 
surface region of the double folding potential can in fact 
be well parametrized by the Woods-Saxon form with the 
diffuseness parameter of around 0.63 fm, and such phe- 
nomenological potential has been as successful as the 
double folding potential. 

In recent years, many experimental evidences have ac- 
cumulated that show that the double folding potential 
fails to account for the fusion cross sections at energies 
close to the Coulomb barrier. That is, the double fold- 
ing potential (and also the Woods-Saxon potential which 
fits elastic scattering) overestimates fusion cross sections 
at energies both above and below the Coulomb barrier, 
having an inconsistent energy dependence to the exper- 
imental fusion excitation function 0, 0, @, H, 0, QU [H| • 
This trend is in accordance with the more recent mea- 
surements of fusion cross sections at extreme subbarier 
energies, that show a much steeper fusion excitation func- 
tions as compared with theoretical predictions [HJ . 

Notice that the scattering process is sensitive mainly 
to the surface region of the nuclear potential, while the 
fusion reaction is also relatively sensitive to the inner 
part. The double folding potential and the Woods- 
Saxon potential are reasonable in the surface region [l3| . 
However, it is not obvious whether they provide reason- 
able parametrizations inside the Coulomb barrier, where 
the colliding nuclei significantly overlap with each other 
[1> EH EH- This is so, particularly because the double 
folding potential takes into account only the so called 
knock-on exchange effect, ignoring all the other exchange 
effects originating from the anti-symmetrization of the 



total wave function of the colliding system jl 61 ] - 

The purpose of this paper is to investigate the radial 
shape of the internucleus potential inside the Coulomb 
barrier and discuss its deviation from the conventional 
parametrizations. To this end, we apply the inversion 
formula based on the WKB approximation [I?) and de- 
termine the internuclear potential directly from the ex- 
perimental data without assuming any parametrization. 
This method was used many years ago by Balantekin et 
al. [l8l |. They assumed a one-dimensional energy in- 
dependent local potential, and found that the inversion 
procedure leads to an unphysical multi-valued potential 
for heavy systems. This analysis has provided a clear 
evidence for inadequacy of the one-dimensional barrier 
passing model for heavy-ion fusion reactions, and has 
triggered to develop the coupled-channels approach. 

Although the analysis of Balantekin et al. is impor- 
tant as it has clarified the dynamics of subbarrier fu- 
sion reactions, it is not satisfactory from the point of 
view of determination of the internucleus potential. Since 
the experimental evidences for inadequacy of the double 
folding potential inside the Coulomb barrier is increas- 
ingly accumulating, it is intriguing to revisit this prob- 
lem by taking into account the current understanding of 
subbarrier fusion. In this connection, we mention that 
the main reason why Balantekin et al. obtained the un- 
physical internucleus potentials is that they did not take 
into account the channel coupling effect, which has by 
now been well understood in terms of barrier distribu- 
tion @, 0, HH HI HI- Our idea here is to apply the 
inversion procedure only to the lowest barrier in the bar- 
rier distribution assuming that the fusion cross sections 
are determined only by it at deep subbarrier energies. We 
will demonstrate below that the internucleus potentials 
thus obtained are well behaved and show a significant 
deviation from the conventional Woods-Saxon shape. 

For a single channel system with a potential V(r), the 
inversion formula relates the thickness of the potential, 
i.e., the distance between the two classical turnin g po ints 
at a given energy E, with the classical action S as [13, EH 

t(E) = r 2 (E)- n (E) (1) 
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where n is the reduced mass between the colliding nuclei 
and Vb is the height of the potential. The classical action 
S(E) is given by 



Jn(E) V n 



(V(r)-E), 



(3) 



and can be obtained once the penetrability P(E) is found 
in some way using the WKB relation P(E) = 1/[1 + 

e 2S(E)^ 

In heavy-ion fusion reactions, it is well known that 
the s-wave penetrability for the Coulomb barrier can be 
approximately obtained from the fusion cross section crf us 

as hi m hi m 



(4) 



where the effective moment of inertia R may depend on 
energy [l8|, This formula assumes that the num- 

ber of classical turning point is two for all partial waves, 
and thus implicitly assumes a deep internuclear poten- 
tial. In the previous application of the inversion formula 
by Balantekin et al, they assumed that the penetrabil- 
ity so obtained was resulted from the penetration of a 
one dimensional energy independent potential (l8j . In- 
stead, here we assume that the penetrability P is given 
as a weighted sum of contribution from many distributed 
barriers, where the distribution arises due to a coupling 
of the relative motion between the colliding nuclei to nu- 
clear intrinsic degrees of freedoms such as collective vi- 
brational or rotational excitations. In this eigen channel 
picture, the penetrability is given by, 



P(E) = J2w n P n (E), 



(5) 



where P n is the penetrability for the n-th eigen-barrier 
and w n is the corresponding weight factor. This concept 
has been well established by now from the experimental 
measurements for the barrier distribution 0, |2l[ as well 
as from numerical calculations of coupled-channels equa- 
tions [H, [24j. In principle, the weight factors w n depend 
on energy if the excitation energy for the intrinsic motion 
is not zero. However, the energy dependence is shown to 
be weak [2-ij . and we assume in this paper that the weight 
factors are energy independent. 

At energies below the lowest eigen barrier (i.e., the 
adiabatic barrier) in the barrier distribution, one expects 
that only the lowest barrier contributes to the total pen- 
etrability, 



P{E) » w P {E). 



(6) 



This indicates that one can apply the inversion formula 
to the lowest eigen barrier using fusion cross sections at 



deep subbarrier energies, after correcting the weight fac- 
tor. In order to demonstrate how this works, Figs. 1(a) 
and 1(b) show the second and the first derivatives of the 
measured Eaf UB Q for the 16 0+ 144 Sm reaction, respec- 
tively. The former quantity is usually referred to as the 
fusion barrier distribution [2l[ . We use the point dif- 
ference formula with AE c _ m =1.8 MeV to carry out the 
derivatives. For this system, one can clearly recognize 
that the barrier distribution has a double peaked struc- 
ture. Correspondingly, the first derivative d(E<jf us ) / dE 
appears to have two steps as a function of energy. If one 
neglects weak coup ling s to the double phonon state in 
the 144 Sm nucleus [25|, the double peaked structure of 
the barrier distribution can be interpreted to originate 
mainly from the coupling to the first 3~ state in 144 Sm. 
Notice that in general a barrier distribution has a peak 
at the energy equal to the height of a potential barrier. 
Assuming that the main peak of the barrier distribution 
around E c , m , ~60 MeV consists only of the contribution 
from the lowest eigen barrier, we scale the first deriva- 
tive d(E(jf us )/dE so that it has a value of 0.5 at the peak 
energy, which we assume to be identical to the position 
of the lowest barrier, The scaling factor corresponds 
to the product of the geometrical factor ttR 2 and the 
weight factor wq (see Eqs. d3| and (J6j) ) . The function 
thus obtained is shown by the filled circles in Fig. 1(c). 
This function can be interpreted as the penetrability for 
the lowest barrier, to which one can apply the inversion 
formula to determine the radial shape. 

The inversion formula yields only the barrier thickness, 
t(E), and one has to supplement either the outer or the 
inner turning points to determine the radial shape of the 
potential [18J. We estimate the outer turning point T2{E) 
using the Coulomb interaction of point charge and the 
Woods-Saxon nuclear potential, 



V N (r) = — 



V 



1 + exp[(r - R )/a] ' 
with the range parameter of 
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and the diffuseness parameter of a=0.63 fm. We adjust 
the depth Vq in order to reproduce the barrier height V& 
determined from the peak position of the barrier distri- 
bution. Since the Coulomb term dominates at the outer 
turning point, except for the region near the barrier top, 
the inverted potential is insensitive to the actual shape of 
nuclear potential employed to estimate the outer turning 
point. The Woods-Saxon potential ([7]) determines not 
only the outer turning point but also the position of the 
potential barrier, Rt,. Following Ref. [Hj], we use 
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R(E) = -{R b + Z P Z T e 2 /E) , 



(9) 



for the effective moment of inertia in Eq. Q . The pene- 
trability obtained by taking into account the energy de- 
pendence of the effective moment of inertia is denoted 



3 




65 
(MeV) 



FIG. 1: Fig. 1(a): The experimental fusion barrier distribu- 
tion for the 16 0+ 144 Sm reaction defined as d 2 (Ea fuB )/dE 2 . 
The experimental data are taken from Ref. JjJ. Fig. 1(b): 
The first derivative of Ea{ UB for the 16 0+ 144 Sm reaction. Fig. 
1(c): The same as Fig.l (b), but normalized so that it is 0.5 at 
the energy of the lower peak in the barrier distribution shown 
in Fig. 1(a) (the filled circle). The open circles are obtained 
in a similar way, but by taking into account the energy de- 
pendence of the effective moment of inertia as given by Eq. 



by the open circles in Fig. 1(c). Although the differ- 
ence between the filled and open circles is not large, 
especially at energies below the barrier, the penetrabil- 
ity behaves slightly better if one considers the energy 
dependence of moment of inertia, since it is saturated 
at unity at high energies. In the actual application of 
the inversion formula shown below, we smooth the data 
points with a fifth-order polynomial fit to the function 
ln[i?crf us /7ri?(£') 2 ] [l8j]. We have confirmed that the re- 
sults do not significantly change even if we use a higher 
order polynomial fit. We also fit the lowest peak of the 
barrier distribution using the Wong formula [26j in order 
to accurately estimate the barrier height V&. 

We have confirmed the accuracy of the inversion proce- 
dure using the theoretical fusion cross sections obtained 
by the computer code CCFULL (27J. For this purpose, we 
consider the 16 0+ 144 Sm system, and generate the fusion 



cross sections by taking into account the excitation to 
the first 3~ state in 144 Sm. We use the same parameters 
as in Ref. [28|. We find that the resultant inverted po- 
tential closely follows the adiabatic potential obtained by 
diagonalising the coupling Hamiltonian at each position 
r [23l [23|. This evidently justifies our procedure for the 
potential inversion discussed above. 




FIG. 2: (Color Online) The adiabatic potential for the 
16 0+ 144 Sm reaction obtained with the inversion method. The 
dashed line is a barrier due to a phenomenological Woods- 
Saxon potential. 



We now invert the experimental data for the 
16 0+ 144 Sm system shown by the open circles in Fig. 1(c) 
in order to obtain the radial dependence of the adiabatic 
barrier. The result of the inversion method is shown in 
Fig. 2. The uncertainty of the inverted potential is esti- 
mated in the same way as in Ref. [Hj]. The dashed line 
shows the barrier due to the Woods-Saxon potential ([7]) 
used to estimate the outer turning points. One clearly 
sees that the inverted potential is much thicker than the 
phenomenological potential at low energies, although it is 
close to the phenomenological potential at energies close 
to the potential barrier. This trend is opposite to what 
Balantekin et al. found in the previous analysis. If there 
was an unresolved peak in the barrier distribution below 
the main peak, one would obtain a much thinner barrier 
than the phenomenological potential, as in the previous 
analysis. We have actually obtained such unphysical thin 
barriers for the 17 + 144 Sm and 16 + 148 Sm systems, 
where the main peak of the barrier distribution is not 
expected to correspond to the lowest eigen barrier [f|. 
Having a thick barrier, rather than a thin barrier, we 
are convinced that the main peak of the barrier distribu- 
tion for the 16 0+ 144 Sm reaction indeed consists of the 
lowest eigen barrier. In this way, the potential inversion 
method could also be used to judge whether there is an 
unresolved barrier in the barrier distribution. 

For the 16 0+ 144 Sm reaction, the experimental data 
exist only from £? cm .=56.6 MeV (that is, about 3.6 MeV 
below the lowest barrier height), and the inverted poten- 
tial has a large uncertainty below this energy. In order 
to discuss the behaviour of the potential far below the 
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barrier, we next consider the 16 O+ 208 Pb reaction, for 
which the fusion cross sections were recently measured 
at deep subbbarier energies [2^|. Figure 3 shows the po- 
tential barrier for this system obtained with the inver- 
sion method. One sees that the behaviour of the poten- 
tial is qualitatively similar to that for the 16 0+ 144 Sm 
reaction shown in Fig. 2. Namely, the inverted poten- 
tial is close to a phenomenological potential in the re- 
gion near the barrier top, but it deviates largely at deep 
subbarrier energies, where the thickness is much larger 
than the phenomenological potential. The thicker the 
potential is, the smaller the penetrability is, and also the 
stronger the energy dependence of the penetrability is. 
The thick potential barrier obtained for the 16 0+ 144 Sm 
and 16 O+ 208 Pb systems is thus consistent with the recent 
experimental observations [13, H§] that the fusion excita- 
tion function is much steeper than theoretical predictions 
at deep subbarrier energies. Although the present anal- 
ysis does not exclude a possibility of a shallow potential 
flol ]. the present study suggests that the origin of the 
steep fall-off phenomenon of fusion cross section can be 
at least partly attributed to the departure of internuclear 
potential from the Woods-Saxon shape. 

For the 16 O+ 208 Pb system shown in Fig. 3, the devia- 
tion of the inverted potential from the phenomenological 
potential starts to occur at around E = 70.4 MeV. It 
is interesting to notice that this energy is very close to 
the potential energy at the contact configuration esti- 
mated with the Krappe-Nix-Sierk potential [13, [H|. In- 
side the touching configuration, the potential represents 



the fission-like adiabatic potential energy surface. The 
effect of such one-body potential has been considered re- 
cently and is shown to account well for the steep fall-off 
phenomena of fusion cross sections [3ll| . The inverted 
potentials which we obtain are thus intimately related to 
the one-body dynamics for deep subbarrier fusion reac- 
tions. 

In summary, we applied the potential inversion 
method, which relates the potential penetrability to the 
thickness of the potential barrier, in order to investigate 
the radial dependence of the internucleus potential for 
heavy-ion fusion reactions. To this end, we assumed that 
the tunneling is well described by the lowest adiabatic 
barrier at deep subbarier energies, and extracted the pen- 
etrability by combining the experimental barrier distri- 
bution and fusion cross sections. We found that the re- 
sultant potential for the 16 0+ 144 Sm and 16 O+ 208 Pb sys- 
tems is much thicker than a barrier obtained with a phe- 
nomenological Woods-Saxon potential. This indicates 
that the steep fall-off phenomenon of fusion cross sec- 
tions recently observed in several systems can be partly 
accounted for in terms of a deviation of internuclear po- 
tential from the Woods-Saxon shape. 

The lowest peak in a barrier distribution is relatively 
well resolved in general for systems involved with a vibra- 
tional target. It would be an interesting future work to 
apply the inversion method systematically to such sys- 
tems and discuss a global potential for heavy-ion colli- 
sions. Another problem is the dynamical effects after 
the touching configuration. In this paper, we exploited 
a barrier distribution picture for subbarrier fusion, but 
assumed an energy independent potential for each dis- 
tributed potential barrier. It would be an interesting 
work to discuss how the thick potential which we ob- 
tained in this paper is related to dynamical effects such 
as the coordinate dependence of reduced mass, and en- 
ergy and angular momentum dissipations, that are other 
promising origins for the steep fall-off phenomenon of fu- 
sion cross sections @, HH, HH • 
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